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Using a Monte Carlo minimization approach we report the global minima for metal clusters 
modelled by the Sutton-Chen family of potentials with N < 80, where N is the number of atoms. 
The resulting structures are discussed in the light of both experimental and theoretical data for 
clusters of the appropriate elements. 



I. INTRODUCTION 

The structure is probably the most fundamental prop- 
erty of a cluster and is important for understanding all as- 
pects of chemical and physical behaviour. Unfortunately, 
there is, as yet, no direct method for the structural de- 
termination of free clusters in molecular beams. Instead, 
one has to measure properties which depend upon ge- 
ometry and then try to infer the structure by compar- 
ing the results with the predictions of models. This ap- 
proach has beenj-combined with techniques such as elec- 
tron diffraption,EJ mass, spectral abundancesp chemical 
reactivityjj magnetisms and x-ray spectroscopy.^ 

For transition metal clusters a wealth of structural in- 
formation is now becoming available from increasingly 
sophisticated experiments. One of the most powerful 
techniques is the flow-reactor approach which probes the 
chemical reactivity of size-selected clusters. For example, 
this method has been applied to nickel clusters using ni- 
trogen as the chemical pwbe to give detailed information 
for all sizes up to A^TIITIj From these data it has been 
possible to make structural assignments arouiid N—\i 
and 55 (sizes of complete Mackay icosahedraEH) , but in 
the size range 29 < N < 48 only one structural assign- 
ment has so far been made because of the large number 
of possible geometries to be considered and the presence 
of multiple isomers.H 

It is, therefore, increasingly important that the theo- 
retician aids the task of experimental interpretation by 
producing reliable structural models. However, for tran- 
sition metal clusters this is an extremely demanding task. 
It is now becoming possible to perform ab initio calcula- 
tions for clusters at the lcucget end of the size range that 
we consider here (N < 80)EJObut only for a few (usually 
high symmetry) geometries. It is not feasible to perform 
the extensive search of the potential energy surface re- 
quired to find the most stable structure. Instead one has 
to use empirical potentials. Even with these simplified 
descriptions of the interatomic interactions it can be a 
very difficult task to find the global minimum for the size 
range considered here. However, we are confident that in 
this work reliable estimates for the global minima have 
been found based upon the -performance of the chosen 
algorithm in previous studies^ and because of the large 



database of structures we have acquired through work on 
clusters bound by the Morse potential.liJIlj 

The three main morphologies that metal and sim- 
ple atomic clusters adopt are icosahedra, decahedra and 
close-packing. Particularly stable examples of each are 
given in Fig. [l]. The icosahedra and decahedra exhibit 
five-fold axes of symmetry which are permitted due to 
the absence— pf translational periodicity. The Mackay 
icosahedror£3 (Fig. |l|b) can be decomposed into twenty 
face-centred cubic (fee) tetrahedra sharing a common 
vertex, and the decahedra are based upon pentagonal 
bipyramids made up of five fee tetrahedra sharing a com- 
mon edge. |T-he most stable decahedral form, the Marks 
decahedrorilj (Fig. 0c), involves further facetting to make 
the cluster more spherical whilst keeping the proportion 
of the surface exhibiting {100} facets to a minimum. The 
most stable fee cluster is the truncated octahedron. All 
three morphologies are commonly seen in metal clusters 
supported on metal surfaces O 




FIG. 1. (a) 38-atom truncated octahedron, (b) 55-atom 
Mackay icosahedron, and (c) 75-atom Marks decahedron. 
These clusters have optimal shapes for the three main types 
of ordered packing seen in this study: face-centred cubic 
(fee), icosahedral and decahedral, respectively. The latter 
two morphologies cannot be extended to the bulk because 
they possess five-fold axes of symmetry. All three structures 
are global minima for the SC 12-6 and 9-6 potentials, and 
only the Mackay icosahedron is not the global minimum for 
the 10-8 potential. 



In this paper we seek to further our understanding of 
transition metal clusters and aid structural assignments 
from experimental data by performing global optimiza- 
tion for clusters interacting via the Sutton-Chen family of 
potentials for N < 80. In section [il] we describe the meth- 
ods used and in section III we describe the structures of 
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the global minima and compare them to experimental 
results and previous theoretical studies. 



II. METHODS 



A. Potentials 



The Sutton-Chen (SC) potential has the form:! 



B = e H 
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c is a dimensionless parameter, e is a parameter with di- 
mensions of energy, a is the lattice constant, and m and 
n are positive integers with n > m. We use_the n, m 
and c parameters given by Sutton and ChenEj for the 
metals Ag, Ni and Au; Rh has the same scaled parame- 
ters as Ag, Cu the same as Ni and Pt is the same as Au, 
so the corresponding results for these metals can simply 
be obtained from their partners by rescaling. For Ag 
and Rh n = 12, m = 6 and c = 144.41, for Ni and Cu 
n = 9, m — 6 and c = 39.432, for Au and Pt n = 10, 
m = 8 and c = 34.408. In the present calculations we 
employed reduced units with e — I and a = 1. The tab- 
ulated results may therefore easily be rescaled for any 
of the above elements. The appropriate energy is given 
simply by multiplying the reduced energy by e whilst the 
coordinates need to be multiplied by a, i.e. the lattice 
constant. The Sutton-Chen potential provides-a reason- 
able description of various bulk properties f3/H3 with an 
approximate many-body representation of the delocal- 
ized metallic bonding. However, it does not include any 
directional terms, which are likely to be important for 
transition metals with partially occupied d shells. 



B. Global Optimization methods 

The main method we used to find the lowest min- 
ima is based- upon Li and Scheraga's Monte Carlo 
minimizationcJ or 'basin-hopping' algorith m .. whi ch we 
have recently explored for several systems .EjOEI This 
approach belongs to the family of 'hypersurface deforma- 
tion' methodsca where the energy is transformed, usually 
to a smoother surface with fewer minima. The lowest 
minimum of the new surface is then mapped back to the 
original surface, but there is no guarantee that the global 
minima on the two surfaces are related and often they 



are not. Ill In contrast, the transformation that we ap- 
ply is guaranteed to preserve the global minimum. The 
transformed energy E is defined by: 



E(X) 



where X represents the vector of nuclear coordinates and 
min signifies that an energy minimization is performed 
starting from X. 

The topography of the transformed surface is that of 
a multi-dimensional staircase. Each step corresponds to 
the basin of attraction surrounding a particular mini- 
mum (the set of geometries where geometry optimiza- 
tion leads to that minimum). The transformation has a 
significant effect on the dynamics. Not only are transi- 
tions to a lower energy minimum barrierless, but they 
can also occur at any point along the boundary between 
basins of attraction, whereas on the untransformed sur- 
face transitions can occur only when the system passes 
along the transition state valley. As a result intrawell 
vibrational motion is removed and the system can hop 
directly between minima at each step. The success of 
the basin-hopping method for potential energy surfaces 
which exhibit multiple funnels has been explained else- 
where in texxns of the thermodynamics of the transformed 
landscaped Shpila*, methods have been used in studies 
of biomolecules.E3~E3 

To explore the E surface we have used canonical Monte 
Carlo (MC) sampling at temperatures of T*=30, 5 and 10 
for the 12-6, 10-8 and 9-6 potentials respectively, where 
the reduced temperature is kT/e. To restrict the config- 
uration space to bound clusters we reset the coordinates 
to those of the current minimum in the Markov chain at 
each step. This objective can also be achieved by placing 
the cluster in a container. 

In our recent application to LJ clusters the MC min- 
imization approach outperformed all other methods in 
the literature, finding all the known lowest energy LJ 
clusters up to 110 atoms including those with non- 
icosahedral global minima. t3 All our published results 
are available in downloadable form from the Cambridge 
Cluster Database.c3 

In a recent study we found the global minima of 
clusters interacting with a Morse potential for all sizes 
up to 80 atoms as a function of the parameter in 
the Morse pptpitial which determines the range of the 
interactions x3tB For this system there are at least 350 
different global minima in this size range. Most of the 
global minima have icosahedral, decahedral or closed- 
packed structures, and these were mainly found by con- 
sidering the structures that have the largest number of 
nearest-neighbour contacts for each morphologyEZI To 
complement the basin-hopping calculations we reopti- 
mized all the Morse global minima and low-lying struc- 
tures for the SC potentials. 
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III. RESULTS 

The energies and point groups of all the global min- 
ima are given in Table ||. The basin-hopping algorithm 
found 95% of the 12-6 global minima, 94% of the 9-6 
and 85% of the 10-8. Five runs from different random 
starting points consisting of 5000 MC steps each were 
used for each cluster size. These results generally confirm 
the utility and robustness of the basin-hopping approach. 
Most of the failures occur when the global minimum has 
a close-packed geometry which is only marginally lower 
in energy than an icosahedral or decahedral structure. In 
these cases the topography of the potential energy sur- 
face is likely to have multiple funnels, a scenario which, 
often makes global optimization extremely difficult E3tB 
We have not made any systematic attempt to find the op- 
timal temperatures or number of steps for the MC runs 
in the present work. 

It is interesting to note that 91% of the 12-6 global 
minima are also global minima for Morse clusters, 63% of 
the 9-6 and 39% of the 10-8. The decreasing percentages 
reflect the relative propensity of these potentials to give 
clusters with ordered structures of the icosahedral, dec- 
ahedral or close-packed morphologies. These values also 
confirm our suggestion that the database of Morse global 
minima would be useful in providing candidate striic-, 
tures for studies with more sophisticated potentials.t3 
Furthermore, the results imply that some of the factors 
determining which isomer of a particular morphology is 
most stable are the same for a simple pair potential and 
for the more realistic many-body potentials used here. 
Indeed virtually all the Sutton-Chen clusters that have 
an ordered icosahedral, decahedral or close-packed struc- 
ture were found in the reoptimization of the larger Morse 
database (global minima and low- lying isomers), where 
the primary structural principle is the maximization of 
the number of nearest neighbours. 

The results presented here match or better all pre- 
viously published results for Sutton-Chen clusters. We 
seem to have found the same structures as Nayak et al. 
who examined all 9-6 clusters with N < 23 (as no en- 
ergies— were given in that paper we cannot be sure of 
this).E3 Uppenbrink and Wales studied a selection of 12-6 
and 10-8 clusters in terms of their thermodynamics and 
dynamics .lJ They attempted to find the global minimum 
by performing regular minimizations from a molecular 
dynamics trajectory. It is interesting to note that for the 
two larger sizes they considered (A=40 and 55) the true 
global minimum was found in only one of the four cases. 



A. Silver and rhodium (SC 12—6) clusters 

The structures of the 12-6 global minima are illus- 
trated in Fig. H and the energies are plotted in Fig. ||a. 
A function of the form a + bN 1/3 + cN 2/3 + dN, which 
has been fitted to the energy, has been subtracted to 



emphasize the size dependence. In Fig. ||b A2-E(iV) = 
E(N+1)+E(N-1)-2E(N) is also illustrated; A 2 E mea- 
sures the stability of a structure with respect to neigh- 
bouring sizes; peaks in A 2 E have been found to correlate 
well with the magic numbers (sizes at which clusters are 
particularly abundant) observed in mass spectra.EiJ 
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FIG. 3. (a) Energies and (b) A 2 E for 
SC 12-6 clusters. In (a) the energy zero, 

E Ag = 940. 6021-994.8068A 1/3 +1126. 5506A 2/3 -1201.39517V , 
where the coefficients have been chosen to give the best fit 
to the energies. 

From Fig. ||a it can be seen that the most stable struc- 
tures occut-at sizes corresponding to complete Mackay 
icosahedraE2l (iV=13 and 55). This result was expected; it 
has been estimated from comparisons of stable sequences 
of clusters that the Mackay icosahedra are lowest in ep. 
ergy for up to a few thousand atoms in this system.L^I 
This result also agrees with an ab initio study of sil- 
ver clusters where, at the few sizes consideced, icosahe- 
dra were always lower than fee structures.!^ The pre- 
dominance of icosahedral morphologies might also be ex- 
pected from the behaviour of the Lennard-Jones potential 
to which the current functional form bears some resem- 
blance. For Lennard-Jones clusters there are only four 
noiMCOsahedral global minima in the range 13 < N < 
80.13-21 However, for the SC 12-6 potential there are only 
30 icosahedral global minima in this size range; at sizes 
between the complete Mackay icosahedra other types of 
structure become lower in energy. Clearly the many-body 
part of the potential is important in determining the most 
stable structure. 
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Initially the growth sequence is similar to that of 
Lennard-Jones clusters. From N=7 to 13 growth oc- 
curs by capping the 7-atom pentagonal bipyramid and 
leads to the 13-atom icosahedron. The one exception 
is the 8-atom cluster for which a deltahedral dodeca- 
hedron is lower in energy than the capped pentagonal 
bipyramid. Growth on a Mackay icosahedron can oc- 
cur in two ways. In the first, the anti-Mackay overlayer, 
the atoms are added in sites which are hexagonal close- 
packed (hep) with respect to the twenty tetrahedra that 
make up the icosahedron; i.e. for the 13-atom icosahe- 
dron atoms are added to the centres of the faces and to 
the vertices. In the second, the Mackay overlayer, the 
atoms are added in sites which are fee with respect to 
the underlying tetrahedra which leads to the next Mackay 
icosahedron. From 7V=14 to 21 growth occurs in the anti- 
Mackay layer, passing through the stable 19-atom dou- 
ble icosahedron (Fig. ||). However, above this size there 
are many types of competing structure and the global 
minimum changes frequently with size. Not until N=51 
are the structures again uniformly icosahedral leading to 
the complete Mackay icosahedron at 55 atoms, However, 
Mackay icosahedra with one and two faces missing do 
give rise to the shallow minima in the energy plot (Fig. 
|a) at iV=46 and 49. 

The two other morphologies exhibited by the global 
minima are decahedral and close-packed. In Table 1 
the close-packed clusters have been divided into those 
that are hep, those that are fee and those that involve 
a mixture of stacking sequences and twin planes. The 
most stable cluster in this intermediate size range is the 
38-atom truncated octahedron (Fig. ||a). The stability 
of this stmcjtare has recently beep recognized in both 
theoreticalEJIla and experimental!!! work. Its shape is 
close to the ideal Wulff polyhedron, and it is the only fee 
structure that is the global minimum for the Lennard- 
Jones potential in this size range. 

For some of the decahedral clusters the (pseudo)- 
fivefold axis is not always obvious; it is obscured by over- 
layers on the {111} faces surrounding the axis at iV=25, 
30, 45 and 48 for this potential (and at iV=45, 48, 58 for 
the 9-6 potential). Also, two global minima (N=22 and 
23) do not belong to any of the ordered morphologies. 
The 23-atom structure is based on two distorted face- 
sharing icosahedra and has been found before for Morse 
clusters;t3 the 22-atom structure is similar. 

Growth on the 55-atom Mackay icosahedron again be- 
gins in the anti-Mackay overlayer; this leads to the weak 
minima in the energy plot of Fig. ||a at N=58 and 61 
which correspond to complete overlayers on one or two 
faces of the icosahedron. However, decahedral structures 
soon become lower in energy. From iV=63 a growth se- 
quence begins which leads to the 75-atom Marks decahe- 



dron. The latter structure's stability is clear from Fig. 
1 

Unfortunately, it is difficult to make any critical as- 
sessment of the performance and reliability of the 12-6 
SC potential because there has been little theoretical or 
experimental work on the structure of silver or rhodium 
clusters in the size range we consider here. 

B. Nickel and copper (SC 9—6) clusters 

The 9-6 global minima are illustrated in Fig. [| and 
the size-dependence of the energies and A2E are given in 
Fig. ||. The latter figure is very similar to that found for 
the 12-6 clusters (Fig. ||): Mackay icosahedra are again 
most prominent, followed by the truncated octahedron 
and the Marks decahedron. However, the differences be- 
tween the depths of the icosahedral (7V=13, 55) and non- 
icosahedral (iV=38, 75) minima in the energy plot in Fig. 

have been reduced compared to Fig ||a; this reflects a 
slight stabilization of the fee and decahedral structures 
with respect to the icosahedra. Moreover, there are no 
subsidiary minima in the energy plot due to icosahedral 
structures at iV=19, 46 and 49. In total there are only 
11 icosahedral minima in the size range 13 < N < 80. 

Two additional close-packed structures at iV=50 and 
59 become more prominent in Fig. ^ The 50-atom clus- 
ter has Dsh point group symmetry and has a twin plane 
passing through the centre; each half of the cluster on 
either side of this twin plane has a structure which is 
a fragment of the 38-atom truncated octahedron. An 
analogous structure is the global minimum for N=79. 
The 59-atom cluster is based on the 31-atom truncated 
tetrahedron with the four sides covered by seven-atom 
hexagonal overlayers. Both are global minima for Morse 
clusters. 

Probably the biggest difference from the 12-6 clusters 
is the large number of global minima in the range N =16 
to 44 which do not belong to any of the ordered morpholo- 
gies. There are 25 such structures. The central atom 
of the 15-atom cluster is 14-coordinate. This structu^a 
is one of the Frank-Kasper coordination polyhedraj2M-3 
and is the. global minimum for a lone=ranged Morse 
potential. Ej A single negative disclinatior£3 runs through 
this cluster. The 16- and 17-atom structures are based 
on the 13-atom Ino decahedron with three and four of 
the square faces capped, respectively. However, the 
structures distort so that two caps approach to form a 
nearest- neighbour contact. They are related to the undis- 
torted decahedta by a single diamond-square-diamond 
rearrangement.^ The structures for iV=21-24 are simi- 
lar to the 23-atom 12-6 global minimum which is made 
up of two distorted face-sharing icosahedra. 
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The 25-atom cluster is based on a C% v fee structure, 
but with the triangular face twisted to remove {100} 
facets. The 34-atom cluster resembles the decahedral 33- 
atom structure; however, the additional atom causes one 
part of the cluster to distort. The 40- atom global mini- 
mum is based on an icosahedral structure with a Mackay 
overlayer (as for 12-6). However, a low coordination 
number atom in that structure has been absorbed into 
the surface. Similarly, the 56-atom and 57-atom clusters 
are based on the Mackay icosahedron. The first adatom 
adsorbs into the surface layer with an accompanying dis- 
tortion, rather than occupying a three-coordinate site on 
the surface. The second adatom then occupies a four- 
coordinate site on the resulting distorted surface. Similar 
avoidance of structures with atoms of low coordination 
number has been_seen in calculations for nickel by Wet- 
zel and DePristoHO The 62- and 63-atom structures are 
based on the 52-atom global minimum, which is a Mackay 
icosahedron with an edge removed. Two triangular faces 
are added over this missing edge. 




FIG. 5. (a) Energies and (b) A 2 E for SC 9-6 clusters. 
In (a) the energy zero, 

E m = 271.8994-292.88737V 1/3 +260.6812iV 2/3 -292.9018iV, 
where the coefficients have been chosen to give the best fit 
to the energies. 

How can these relatively disordered structures become 
global minima, and why do they disappear at larger sizes? 
Probably because most of these geometries are more 
spherical, have a larger number of nearest-neighbour con- 
tacts and involve fewer low coordination number atoms 
than the possible structures based on one of the or- 



dered morphologies. However, this advantage is coun- 
terbalanced by the introduction of considerable strain. 
As the eaeijgetic penalty for the strain scales with the 
volume,OE3 it is only possible to accommodate it at 
smaller sizes. 

For the elements described by this potential, there is a 
wealth of theoretical and experimental work to compare 
with our results. This fortunate situation is mainly due 
to the pioneering work of Riley and coworkers on nickel 
clusters. First we compare our results to other theoreti- 
cal studies which use different descriptions of the interac- 
tions. Of these studies the most comprehensive is that of 
DePristo and coworkers who used the corrected effective, 
medium theory to look at all clusters up to A=55E3't3 
They found that at iV=13 and 55 the favoured geome- 
tries were icosahedral. However, for many of the interme- 
diate sizes, especially when the icosahedral structure at 
that size would involve a low coordination number sur- 
face atom, the global minima did not belong to any of 
the ordered morphologies. At some sizes these disordered 
structures appear to be the same as found in this study, 
e.g. iV=15, 18, and 25. DePristo et al. did not identify 
any of the global minima as close-packed or decahedral, 
although it seems that their 38-atom structure is in fact 
a distorted truncated octahedronS 

Using a tight-binding model Lathiotakis et al. rStud- 
icd a selection of clusters with up to 55 atoms .c2l Of 
the structures they considered, they found that Mackay 
icosahedra were lowest in energy at A=13 and 55, but 
at some intermediate sizes fee structures were lower in 
energy than icosahedral clusters. However, only a few 
structures were considered at each size and these may 
not include the global minima. 

Montcjano et al. used an embedded-atom potential to 
compare the energies of icosahedral structures, in partic- 
ular to locate the size at which the transition from an 
anti-Mackay to a Mackay overlayer occurredo However, 
these results are not relevant to many sizes since the- 
ory and experiment suggest that other non-icosahedral 
structures are lowest in energy. 

Finally, Hu et al. attempted to model the interatomic 
interactions using a long-ranged Morse potential. e2I As. 
expected for a Morse potential with range parameterEj 
po=3.54 they observed icosahedral structures with an 
anti-Mackay overlayer up to 7V=40. However, a Morse 
potential is probably not a good approximation to the 
real nickel potential; therefore these results differ from 
experimental and other theoretical studies. 

Icosahedral structure was first identified for nickel clus- 
ters by looking at the size dependence off-the chemical 
reactivity with various probe molecules. cJlj Features 
were found at sizes corresponding to complete Mackay 
icosahedra, and also for icosahedral structures with sta- 
ble surface overlayers. Subsequently, icosahedral magic 
numbers [have also been seen in mass spectra for 50 < 
N < 800£3 

The most detailed information on the structure at par- 
ticular sizes has come from work using nitrogen probe 



7 



molecules B~a This technique has enabled structural as- 
signments to be made for the majority of clusters in the 
range N < 28, 49 < N < 71 and at A=38. However, in 
the intermediate size range, structural assignments based 
on these nitrogen experiments have not yet been pub- 
lished due to the difficulty in interpreting the experimen- 
tal data. 

In agreement with our work Parks et al. found that 
the structures at A=13 and 55 are Mackay icosahedra, 
whilst iV=38 is a truncated octahedron. In the range 
13 < A < 26 the experiments with nitrogen indicate 
that most of the structures are formed by the addition 
of an anti-Mackay overlayer to the 13-atom icosahedron, 
from N—A9 to 55 the Mackay overlayer is completed, 
and then up to N=71 an anti-Mackay overlayer grows 
on the 55-atom Mackay icosahedron.,— The reactivity of 
nickel clusters with ammonia, water, til hydrogerJla and 
deuteriumcH also suggests that icosahedral structure per- 
sists above N=7l. These assignments differ from our 9-6 
results in that the Sutton-Chen potential seems to under- 
estimate the stability of icosahedral structures. However, 
non-icosahedral clusters do seem to occur in the range 
29 < N < 48, in agreement with the SC 9-6 potential. 
If disordered global minima are indeed observed in this 
size range for real nickel clusters this would help to ex- 
plain the difficulties experienced by Riley and coworkers 
in assigning structures. 

To further facilitate comparisons with experiment we 
have estimated the number of binding sites for nitrogen 
using the rules formulated by Parks et al. These are: (1) 
N2 binds directly to the nickel atoms; (2) a nickel atom 
with a coordination number of four or less binds two N2 
molecules; (3) nickel atoms with a coordination number 
of five to eight will readily bind one N2 molecule; (4) 
nickel atoms with a coordination number of nine bind 
N2 molecules weakly or not at all; and (5) nickel atoms 
with a coordination number of ten or more do not bind 
N2 . To clarify rule (4) it has been found that for smaller 
clusters N2 can bind to a nine-coordinate atom, but for 
the larger clusters (N > 49) no evidence has been found 
for this type of binding. The number of binding sites 
and the coordination number analysis for the 9-6 global 
minima are given in Table [n]. For the clusters with an 
ordered morphology the differentiation between nearest 
neighbours and next-nearest neighbours is clear. How- 
ever, this differentiation becomes ambiguous for the dis- 
ordered structures and so the coordination number anal- 
ysis, and sometimes the number of N2 binding sites, be- 
comes dependent on the nearest-neighbour criterion that 
is used. 

The magnetic moments of size-selected nickel clustery 
have been measured for all clusters with up to A=200.Q 
For small N there is considerable variation of the mag- 
netic moment with size. It was concluded that features at 



A=13 and 55 indicate icosahedral structured However, 
it is difficult to decipher the structural information that 
is contained in other features of the magnetic moment 
size-dependence. 

Much less structural information is available for cop- 
per clusters. Recent experiments suggest that icosahedra 
predominate up to N ~ 25QO. atoms and above this fee 
clusters are more prevalento Although that study en- 
compasses clusters far larger than we consider here, our 
results are not inconsistent with this finding. 

C. Gold and platinum (SC 10—8) clusters 

The 10-8 global minima are illustrated in Fig. ^ and 
the size-dependence of the energies and A2E are given in 
Fig. @. The latter fi gure shows a very different pattern 
from that seen for the 12-6 or 9-6 clusters: there are 
no signatures due to icosahedral clusters. Although the 
13-atom icosahedron is the global minimum, above this 
size there are no global minima with ordered icosahe- 
dral structures. Instead Fig. ^a is dominated by features 
due to particularly stable close-packed and decahedral 
clusters. The 38-atom fee truncated octahedron and the 
75-atom Marks decahedron exhibit the deepest minima 
111 Fig. 0a. There are also minima in the energy plot at 
7V=64 and 71 due to incomplete Marks decahedra, and 
at A=50, 61 and 79 due to close-packed clusters. The 
50-atom global minima is the 'twinned truncated octa- 
hedron' that is the global minimum for the 12-6 an 9-6 
potentials; the 61-atom structure is a fragment of the 
79-atom twinned truncated octahedron that is the global 
minimum for the 12-6 and 9-6 potentials, and the 79- 
atom global minimum is a truncated octahedron. That 
the 79-atom global minimum is fee illustrates the 10- 
8 potential's greater preference for fee structures rather 
than close-packed structures involving twin planes. For 
A=54-56 fee structures are also lowest in energy and 
the 59-atom close-packed structure with Td symmetry is 
not the global minimum in contrast to the 12-6 and 9-6 
potentials. 

As for the 9-6 potential, between A=13 and 55 there 
is a tendency to form structures that do not fit neatly 
into one of the ordered morphologies. In total there are 
31 such structures compared to 25 for the 9-6 poten- 
tial. The 14- and 21-atom structures are based on the 
SC 9-6 15-atom global minimum. The 14-atom cluster 
has one of the 7-coordinate vertices removed and the 21- 
atom structure is formed when 6 atoms are added to the 
faces surrounding, a vertex, thus extending the negative 
disclination line.ll3 The structures for A=15-20, like the 
A=15, 16 SC 9-6 global minima, are based on distorted 
decahedra. Similar distorted decahedra are also observed 
at A=35, 37 and 41. 
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FIG. 7. (a) Energies and (b) A 2 E for 

SC 10-8 clusters. In (a) the energy zero, 

Eau = 277.5753-276.0472iV 1/3 +192.0783iV 2/3 -305.9338iV, 
where the coefficients have been chosen to give the best fit 
to the energies. 



The global minima for iV=22-30 all structurally re- 
lated. The smaller clusters resemble the 23-atom D%h 
SC 12-6 global minimum and the larger ones are related 
to a 30-atom structure which is made up of three inter- 
penetrating Dzh units. Interestingly, the 24- and 30-atom 
structures produce minima in the energy plot in Fig. 0a. 

Many of the other disordered structures have, in part, 
the surface structure of an incomplete Mackay icosahe- 
dron, but are distorted in various ways, e.g. N=31, 36, 
42, 44, 46, 48, 52 and 53. The 53-atom structure is based 
on the Mackay icosahedron with three adjacent vertices 
removed, an atom added in the centre of this incomplete 
face and accompanying distortions. The 52-atom struc- 
ture is then formed simply by the removal of one atom. 

The main work to which our results can be compared 
is that of Whetten and (psworkers on gold clusters pas- 
sivated by alkylthiolates.EjTjj In these experiments the 
clusters selectively form specific sizes, which were isolated 
by fractional crystallization. Each fraction has a narrow 
size distribution. The clusters presumably form these 
specific sizes because of their stability. Furthermore, it is 
thought that the passivating surface layer does not per- 
turb the structure significantly. If this is true then the 
observed structures reflect those of the free clusters. 

Most of the cluster sizes formecL-sKere close to those 
expected for truncated octahedra.EJO For the small- 
est sizes detailed structural investigations were made by 



comparing experimental x-ray diffraction patterns with 
those calculated from structural models. This compari- 
son led to the identification of the fijatlipns with N ~ 75, 
101 and 146 as Marks decahedra,EJ£3 and those-,with 
N ~ 225, 459 as twinned truncated octahedra.Eil Re- 
cently, a fraction with N ~ 38 has also been isolated; the 
diffraction patterns suggest that it is the fee truncated 
octahedronH 

Our results for SC 10-8 agree very well with these ex- 
perimental data. The two clusters observed experimen- 
tally for N < 80 are the two that we find to be the most 
stable, namely the 38-atom truncated octahedron and 
the 75-atom Marks decahedron. 

Cleveland et al. have performed a number of theoret- 
ical studies to interpret the experimental results on pas- 
sivated gold clusters .EilH They used an embedded-atom 
potential and looked at several sequences of structures. 
As in the present study, they found that truncated oc- 
tahedra and Marks decahedra were most stable. They 
also found that introducing a twin plane into a trun- 
cated octahedron leads to an increase in energy, e.g. for 
N=79 the fee Oh structure was slightly lower in energy 
than the twinned D 3 h structure. This trend seems to 
disagree with the experimental results — at N ~ 225, 459 
the diffraction patterns seem to suggest that the struc- 
ture has a twin plane. Garzon and Posado-Amirillas also 
used the,embedded-atom method to look at a 55-atom 
cluster.E3 They found that a disordered structure was 
lower in energy than the icosahedron. 

We also note that at 7V=55 the fee cuboctahedron is 
not the global minimum; it actually lies 78. le above the 
fee global minimum. In connection with this comment, 
it is interesting that a recent reinvestigation of ligated 
55-atom gold clusters, which were originally thought to 
be cuboctaJasdral,E3'E3 seems to disprove this structural 
assignment £9 

The only calculations for platinum that are relevant to 
our results are by Sachdev et aZ.Eafl9 For clusters with 
up to 60 atoms they found that the lowest energy struc- 
tures were disordered. In particular at iV=55 these struc- 
tures were lower than the icosahedron and cuboctahe- 
dron. These results have some overlap with ours but it is 
not clear whether Sachdev et al. found the global minima 
for the larger sizes. 



IV. CONCLUSIONS 

In this paper we have found the likely global minima 
for clusters with up to 80 atoms whose interactions are 
described by the SC family of potentials. These poten- 
tials can be used to model silver, rhodium, nickel, cop- 
per, gold and platinum clusters. The results are encour- 
aging, thus confirming the utility of these empirical po- 
tentials. For example, the most stable structures gen- 
erally appear to agree with experiment. For nickel clus- 
ters the 13-atom and 55-atom Mackay icosahedra, the 38- 
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atom truncated octahedron and the 75-atom Marks dec- 
ahedron are particularly stable. The former three clus- 
ters have been unambiguously identified in experimental 
studies of chemical reactivity using nitrogen as a probe 
molecule, and it would be interesting to see if extensions 
of these experiments to N=75 could identify the Marks 
decahedron. For gold clusters we find the 38-atom trun- 
cated octahedron and the 75-atom Marks decahedron to 
be particularly stable; Whetten and coworkers have been 
able to isolate-both these clusters when passivated by 
alkylthiolates .L3~Ej 

Our results should also be useful in aiding structural 
assignments from experimental data on size-selected clus- 
ters. However, it would be surprising if these empirical 
potentials could accurately reproduce all the intricacies 
of a particular cluster growth sequence. This is espe- 
cially true when the energy differences between compet- 
ing structures are small. In such cases the character of 
the global minimum will be particularly sensitive to the 
accuracy of the potential. 

When comparing the structures described here to ex- 
periment, aside from possible inaccuracies in the poten- 
tial, it should also be remembered that the global mini- 
mum is only rigorously the free energy global minimum 
at absolute zero. At finite temperature entropic effects 
may play a role in determining the most stable struc- 
ture. These entropic effects are most likely to be in- 
fluential when the energy gap between the global min- 
imum and other low energy minima is small. For ex- 
ample, it has been shown that for a 38-atom Lennard- 
Jones cluster the structure changes Jipm fee to icosa- 
hedral as the temperature increases.Ea Similarly, for a 
75-atom Morse cluster with a medium-ranged potential 
there is a transition from a Marks decahedron to icosa- 
hedral structures Ji-3 Both these transitions stem from the 
larger entropy of the icosahedra; there are many icosa- 
hedral minima which have energies just above the global 
minimum, whereas there is a large energy gap between 
the global minimum and the next lowest energy mini- 
mum with the same morphology. However, probably a 
more common effect of temperature when energy differ- 
ences between low energy structures are small will be the 
presence of multiple isomers. 

Finally, the results presented here further illustrate the 
power of the 'basin-hopping' or Mcpte Carlo minimiza- 
tion global optimization algorithmic This method has 
enabled us to locate global minima with some confidence 
for systems with up to 240 degrees of freedom, some of 
which exhibit a multiple funnel potential energy surface 
topography. 
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TABLE I. Global minima for Sutton-Chen 12-6, 9-6 and 10-8 clusters with N < 80. For each minimum the energy and point 
group are given and a structural assignment made if possible. The structural categories are: icosahedral with an anti-Mackay 
(aM) or a Mackay over layer (M); decahedral with n atoms along the decahedral axis (dn); close-packed fee (f), hep (h), or 
involving a mixture of stacking sequences and twin planes (c); and clusters involving disclination lines (dis). 



12-6 9-6 10-8 



N 


Energy/e 


PG 




Energy / e 


PG 




Energy/e 


PG 




3 


-1704.6905 


D 3h 




-480.8560 


D 3h 




-633.7771 


D 3h 




4 


-2601.8447 


T d 




-709.5396 


T d 




-904.1153 


T d 




5 


-3461.3452 


D 3h 




-929.7341 


D 3h 




-1163.7670 


D 3h 




6 


-4378.8875 


o h 




-1163.9640 


o h 




-1433.8252 


o h 




7 


-5271.2947 






-1388.5116 


D 5h 




-1695.8893 


D 5h 




8 


-6129.7564 


D 2d 




-1611.8509 


D 2d 




-1954.1206 


Di d 




9 


-7048.7552 


Civ 




-1839.9790 


Civ 




-2218.1861 


Civ 




10 


-7972.0971 


c 3v 




-2072.2436 


C 3v 




-2481.7019 


c 3v 




11 


-8889.9627 


Civ 




-2302.4939 


Civ 




-2741.5489 


Civ 




12 


-9871.2458 






-2543.1611 


Cr,v 




-3005.9274 


Ci 




13 


-10968.5082 


Ih 


M 


-2808.5765 


Ih 


M 


-3280.3843 


Ih 


M 


14 


-11798.8479 


c 3v 


aM 


-3023.9716 


C 3v 


aM 


-3549.4023 


Clj v 


dis 


15 


-12742.9841 


Civ 


aM 


-3267.5300 


D 6d 


dis 


-3825.6495 


Ci v 




16 


-13672.6475 


Cs 


aM 


-3505.2600 


Cs 




-4101.6928 


Cs 




17 


-14606.3231 


c 2 


aM 


-3742.6166 


Civ 




-4371.8696 


Cs 




18 


-15535.3810 


Cs 


aM 


-3978.1268 


Civ 




-4640.1812 


Cs 




19 


-16595.0561 


D 5h 


aM 


-4221.3539 


D 5h 


aM 


-4906.0562 


C s 




20 


-17510.9209 


Civ 


aM 


-4455.5507 


Civ 


aM 


-5171.2334 


Civ 




21 


-18433.0300 


Ci 


aM 


-4700.7823 


Ci 




-5439.9207 


Cqv 


dis 


22 


-19422.7209 


c a 




-4949.3235 


C s 




-5716.6689 


C s 




23 


-20383.3977 


D 3h 




-5191.4317 


c 2 




-5990.7388 


Ci 




24 


-21315.4208 


Civ 


h 


-5427.4229 


Ci 




-6266.6622 


Cs 




25 


-22339.6319 


C 3v 


d3 


-5670.7723 


C 3 




-6533.1606 


c\ 




26 


-23337.2211 


D 3h 


h 


-5920.3488 


D 3h 


h 


-6803.6959 


Cs 




27 


-24284.3891 


C a 


h 


-6165.3671 


Cs 




-7080.7248 


C s 




28 


-25276.9501 


C 3v 


M 


-6411.2387 


Ci 




-7354.7939 


Ci 




29 


-26263.2779 


Civ 


d4 


-6654.0358 


Ci 




-7632.7382 


Ci 




30 


-27253.8536 


Civ 


d3 


-6903.2657 


Cs 




-7909.1542 


C 3v 




31 


-28274.4371 


Civ 


d4 


-7153.4410 


c 3 




-8171.0816 


Ci 




32 


-29265.3320 


Civ 


M 


-7400.8234 


Did 




-8451.6848 


c 3 




33 


-30274.9603 


Civ 


d4 


-7644.6441 


Cs 


d4 


-8726.4506 


Cs 


d4 


34 


-31231.7697 


Civ 


c 


-7889.3674 


c 2 




-9005.3502 


c 3 




35 


-32280.3945 


Civ 


d4 


-8147.0475 


D 3 




-9276.2927 


Cs 




36 


-33253.9352 


Cs 


M 


-8392.4962 


Civ 




-9551.6256 


Civ 




37 


-34302.6067 


C 3v 


c 


-8646.8835 


C 3v 


c 


-9836.6867 


Civ 




38 


-35419.9804 


o h 


f 


-8917.7056 


o h 


f 


-10117.2454 


o h 


f 


39 


-36364.8587 


C&v 


f 


-9156.5715 


C&v 


f 


-10389.6477 


C^v 


f 


40 


-37324.3708 


Cs 


M 


-9397.0850 


Cs 




-10661.9303 


D 4h 


f 


41 


-38316.5698 


Cs 


c 


-9643.6606 


Cs 


c 


-10934.9766 


Cs 




42 


-39301.6696 


C s 


M 


-9892.1913 


C s 




-11211.5400 


C s 




43 


-40341.8543 


Cs 


M 


-10140.5484 


Cs 




-11490.3063 


Civ 


d4 


44 


-41310.9157 


Ci 


M 


-10391.3783 


c 2 




-11767.0685 


Cs 




45 


-42345.0912 


C 3 


d4 


-10642.7040 


C s 


d4 


-12039.6977 


Ci 




46 


-43436.2827 


Civ 


M 


-10900.4123 


Civ 


M 


-12318.3028 


c 3 




47 


-44405.1884 


Ci 


M 


-11145.6538 


Ci 


M 


-12595.3291 


Civ 


d4 


48 


-45470.1069 


Civ 


d4 


-11411.4049 


Civ 


d4 


-12875.3949 


Ci 




49 


-46521.2131 


C 3v 


M 


-11662.0840 


C 3v 


M 


-13150.8235 


D 5h 


d4 


50 


-47518.6719 


D 3h 


c 


-11920.7434 


D 3h 


c 


-13433.4182 


D 3h 


c 


51 


-48522.4267 


Civ 


M 


-12160.5960 


Cs 


c 


-13705.6730 


Cs 


c 


52 


-49616.1377 


C 3v 


M 


-12424.6351 


Civ 


M 


-13986.6288 


Civ 




53 


-50706.4665 


Civ 


M 


-12689.2438 


Civ 


M 


-14261.3449 


C 3v 




54 


-51796.0777 


Cr,v 


M 


-12953.6990 


Cr,v 


M 


-14542.6311 


Civ 


f 


55 


-52884.6806 


I h 


M 


-13217.8963 


h 


i 


-14814.5225 


Ci 


f 
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TABLE I. continued. 



12-6 9-6 10-8 



N 


Energy/e 


PG 




Energy/e 


PG 




Energy/e 


PG 




56 


-53756.6516 




M 


-13445.8961 


Cs 




-15098.5064 


D 2h 


f 


57 


-54700.1733 


c s 


aM 


-13684.1489 


C s 




-15377.3233 


C2v 


d5 


58 


-55753.8515 


Civ 


aM 


-13943.4015 


D 3h 


d4 


-15647.9269 


C s 


d5 


59 


-56751.4572 


T d 


c 


-14202.4032 


T d 


c 


-15928.8952 


Civ 


d5 


60 


-57763.6760 


C a 


aM 


-14445.8412 


C2v 


d5 


-16206.2020 


Civ 


d5 


61 


-58809.0448 


Cl v 


aM 


-14703.3954 


Cz v 


c 


-16495.9638 


Csv 


c 


62 


-59765.2180 


C2v 


aM 


-14951.2644 


Cs 




-16768.1774 


Cs 


c 


63 


-60822.3826 


C a 


d5 


-15213.6371 


C2v 




-17048.9913 


Cs 


c 


64 


-61925.6244 


C2v 


d5 


-15480.8530 


C'2v 


d5 


-17336.8555 


Civ 


d5 


65 


-62903.7387 


Civ 


d5 


-15724.8185 


C2v 


d5 


-17610.0214 


Civ 


d5 


66 


-63959.3105 


Cs 


d5 


-15985.0232 


Cs 


d5 


-17890.9860 


Cs 


d5 


67 


-65011.2767 


C2v 


d5 


-16244.0284 


C2v 


d5 


-18169.0668 


Civ 


d5 


68 


-65980.5983 


Cz v 


c 


-16490.5019 


Cz v 


c 


-18448.5177 


Csv 


c 


69 


-67020.4042 


Ci 


d5 


-16744.2741 


C2v 


d5 


-18725.1157 


Ci 


d5 


70 


-68114.9462 


Cs 


d5 


-17012.0775 


Cs 


d5 


-19008.8021 


Cs 


d5 


71 


-69216.6518 


C2v 


d5 


-17279.8708 


C2v 


d5 


-19300.7255 


Civ 


d5 


72 


-70171.4663 


Ci 


d5 


-17520.4788 


Ci 


d5 


-19573.7638 


Cs 


d5 


73 


-71225.8547 


C2v 


d5 


-17779.2971 


Cs 


d5 


-19855.7668 


Cs 


d5 


74 


-72318.7243 


Cr, v 


d5 


-18047.0929 




d5 


-20138.8921 


C\,v 


d5 


75 


-73421.0521 


D 5h 


d5 


-18315.1577 


D sh 


d5 


-20431.3452 


D 5h 


d5 


76 


-74375.6975 


Cs 


d5 


-18554.6196 


C2v 


d5 


-20704.2050 


Civ 


d5 


77 


-75430.9852 


C2v 


d5 


-18814.0659 


C2v 


d5 


-20978.9269 


Civ 


d5 


78 


-76385.4318 


Ci 


d5 


-19054.3977 


Cs 


d5 


-21252.2801 


Cs 


d5 


79 


-77456.0255 


D 3h 


c 


-19321.0094 


D 3h 


c 


-21547.5979 


O h 


f 


80 


-78414.6271 


Cs 


c 


-19560.8966 


Cs 


c 


-21819.9952 




f 



14 



TABLE II. Coordination number analysis and estimated number of N 2 binding sites for SC 9-6 clusters. The number of 
binding sites has been calculated using the rules given in the text. The values in brackets are appropriate if nine-coordinate 
binding sites are included. For 38 atoms and above there is no evidence of this type of binding and the alternative values have 
been omitted. The nearest-neighbour criterion used is 0.8a. 





N2 binding 








coordination number 






N 


sites 


3 


4 


5 


6 7 8 


9 


> 10 



4 8 4 

5 10 2 3 

6 12 6 



7 


12 


5 


2 












8 


12 


4 


4 












9 


13 


4 


2 


2 




1 






10 


12(13) 


3 


3 


3 






1 




11 


12 


2 


4 


4 








1 


12 


11 




5 


6 








1 


13 


12 






12 








1 


14 


14 






9 


3 






1 


15 


14 






12 


2 






1 


16 


16 


1 


2 


12 








1 


17 


18 


2 




10 


4 






1 


18 


19 


2 




8 


7 






1 


19 


17 






12 




5 




2 


20 


19 


1 




10 


2 


5 




2 


21 


19 




1 


12 


1 


5 




2 


22 


20 




3 


11 


4 


2 




2 


23 


21 






14 


6 


1 




2 


24 


22 




2 


10 


6 


4 




2 


25 


21(22) 






12 


3 


6 


1 


3 


26 


21(24) 




6 


12 




3 


2 


3 


27 


24 






13 


4 


7 




3 


28 


25 






13 


7 


5 




3 


29 


26 




2 


12 


6 


6 




3 


30 


26 




1 


13 


4 


8 




4 


31 


24(27) 




3 


12 


3 


6 


3 


4 


32 


24(28) 






12 






4 


4 


33 


26(26) 




4 


14 


2 


6 


2 


5 


34 


26(30) 




2 


14 


6 


4 


4 


4 


35 


30 






12 


6 


8 




5 


36 


30 






16 


6 


8 




6 


37 


27(31) 




3 


18 


3 


3 


4 


6 


38 


24 






24 






8 


6 


39 


26 


1 




20 


4 




8 


6 


40 


33 






14 


6 


13 




7 


41 


29 




3 


18 


4 


4 


5 


7 


42 


35 




2 


18 


1 


14 




7 


43 


29 




3 


18 


4 


4 


7 


7 


44 


35 






14 


10 


11 




9 


45 


32 




2 


19 


5 


6 


3 


10 


46 


35 






16 


2 


17 


2 


9 


47 


36 




1 


16 


2 


17 


2 


9 


48 


35 






22 


2 


11 


2 


11 


49 


36 






15 


3 


18 


3 


10 


50 


30 






24 


6 




8 


12 


51 


32 


1 




22 


7 


1 


8 


12 


52 


39 






14 


4 


21 




13 


53 


40 






10 


10 


20 




13 


54 


41 






11 


5 


25 




13 


55 


42 






12 




30 




13 



15 



TABLE II. continued. 



N2 binding coordination number 



N 


sites 


3 4 


5 


6 


7 


8 


9 


> 10 


56 


43 






10 


11 


22 




13 


57 


45 


1 




9 


9 


25 




13 


58 


39 






24 


6 


9 




19 


59 


36 






24 


12 




4 


19 


60 


41 






22 




13 


4 


15 


61 


36 






24 


6 


6 


12 


13 


62 


46 






15 


3 


28 




16 


63 


46 






16 




32 




17 


64 


40 






22 


8 


10 


6 


18 


65 


42 


1 




22 


6 


12 


6 


18 


66 


42 




2 


20 


7 


13 


6 


18 


67 


44 






22 


6 


16 


6 


17 


68 


39 






24 


9 


6 


13 


16 


69 


41 






20 


18 


3 


8 


18 


70 


42 






21 


13 


8 


8 


19 


71 


43 






22 


8 


13 


8 


20 


72 


43 


1 




21 


8 


12 


10 


20 


73 


41 




2 


20 


9 


10 


11 


21 


74 


41 






21 


15 


5 


10 


23 


75 


42 






22 


10 


10 


10 


23 


76 


44 


1 




18 


14 


10 


10 


23 


77 


44 




2 


18 


12 


12 


10 


23 


78 


46 


1 




20 


10 


14 


10 


23 


79 


39 






24 


12 


3 


18 


22 


80 


41 


1 




22 


13 


4 


18 


22 



16 



